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Abstract Various methods of helioseismology are used to study the subsurface prop- 
erties of the sunspot in NOAA Active Region 9787. This sunspot was chosen because 
it is axisymmetric, shows little evolution during 20-28 January 2002, and was observed 
continuously by the MDI/SOHO instrument. AR 9787 is visible on helioseismic maps 
of the farside of the Sun from 15 January, i.e. days before it crossed the East limb. 

Oscillations have reduced amplitudes in the sunspot at all frequencies, whereas a 
region of enhanced acoustic power above 5.5 mHz (above the quiet-Sun acoustic cutoff) 
is seen outside the sunspot and the plage region. This enhanced acoustic power has 
been suggested to be caused by the conversion of acoustic waves into magneto-acoustic 
waves that are refracted back into the interior and re-emerge as acoustic waves in 
the quiet Sun. Observations show that the sunspot absorbs a significant fraction of 
the incoming p and / modes around 3 mHz. A numerical simulation of MHD wave 
propagation through a simple model of AR 9787 confirmed that wave absorption is 
likely to be due to the partial conversion of incoming waves into magneto-acoustic 
waves that propagate down the sunspot. 
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Wave travel times and mode frequencies are affected by the sunspot. In most cases, 
wave packets that propagate through the sunspot have reduced travel times. At short 
travel distances, however, the sign of the travel-time shifts appears to depend sensitively 
on how the data are processed and, in particular, on filtering in frequency- wavenumber 
space. We carry out two linear inversions for wave speed: one using travel- times and 
phase-speed filters and the other one using mode frequencies from ring analysis. These 
two inversions give subsurface wave-speed profiles with opposite signs and different 
amplitudes. 

The travel-time measurements also imply different subsurface flow patterns in the 
surface layer depending on the filtering procedure that is used. Current sensitivity 
kernels are unable to reconcile these measurements, perhaps because they rely on im- 
perfect models of the power spectrum of solar oscillations. We present a linear inversion 
for flows of ridge- filtered travel times. This inversion shows a horizontal outflow in the 
upper 4 Mm that is consistent with the moat flow deduced from the surface motion of 
moving magnetic features. 

From this study of AR 9787, we conclude that we are currently unable to provide 
a unified description of the subsurface structure and dynamics of the sunspot. 

Keywords Sun • Sunspots • Helioseismology 



1 Introduction 

One of the main goals of solar physics is to understand the physical processes respon- 
sible for solar magnetism and activity. This requires the study of magnetic flux tubes, 
their transport and dynamics in the convection zone, and their emergence at the solar 
surface in the form of sunspots and active regions. The overall nature of sunspots is 
still a matter of debate. Many open questions remain concerning their structure and, 
above all, their formation and stability. How can regions of such intense magnetic flux 
come into existence and remain stable over several days, weeks, and sometimes months? 
Another common question is whether sunspots are monolithic magnetic fl ux tubes or 
have a spaghetti- like structure ([Parker 1979). Schiiss ler and Rempell (2005) proposed a 
scenario whereby a sunspot expands rapidly below the surface during the early stages 
of its formation, leading to a disconnection from its magnetic roots. This disconnection 
may allow a transition to a spaghetti-like subsurface structure. As to the sta bility of 
sunsp ots, it may be due to the presence of surface and subsurface collar flows ( Parkerl 



ll979T ). Other questions concern the energetics of sunspots, the flow of heat through 
and around sunspots, and the nature of magnet oconve ct ion at kilogauss fields. W hat 
is known about sunspots has been summarized by, e.g., Thomas and Weisd (|l992h and 
ISolankil (|2003l) . For a descr iption of various mag netostatic sunspot models, we refer the 
reader to ljahnl (|l992h and lRempel et al.1 (|2008h . 



In this paper we will discuss the potential of helioseismology to probe the subsurface 
structure of sunspots, with the hope of answering, one day, some of the questions 
listed above. Local helioseismolog y includes severa l met hods of analysis, which have 
been described in some detail by iGizon and Birchl (2005). All these methods rely on 
continuous time series of Doppler images of the Sun's surface. Fourier- Hankel analysis 
was deve loped to study the relationship between ingoing and outgoing waves around a 
sunspot ([Braun et al.lll987l : lBraun| [l995). Ring-diagram analysis consists of analysing 
the frequencies of solar acoustic waves over small patches of the solar surface ([Hillll 1988k 
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lAntia and BasulfeoOTh . Time-distance helioseismol ogy (TP) measures th e travel times 
of wave packets moving through the solar interior (|Duvall Jr. et al.l ll993). Helioseismic 
holography (HH) u ses the observed wave fie ld at the solar surface to infer the wave field 
at different dept hs jLindsev and Braunlll997|). A sum mary of recent results is provided 
bv lGizonl (fcOOlSh and lThompson and Zharkovl (|2008h . 

There ha ve b een sev e ral st udies of sunspots using helioseismology. iBraun et al.l 
<|1987L Il992ah and Braun (1995) used Fourier- Hankel decomposition to measure wave 
absorption and scattering phase shifts caused by sunspots. The absorption is believed 
to be the result o f a partial conversion of incoming p modes into slow magnetoa- 
coustic waves (e.g., ISpruit and Bogda n1 ll992l : ICallv1l2000h . Observ ational signat ures of 
the m ode conversion process have been discusse d, for ex ample, by Sc hunker and Callvl 
(|2006h . Agreement bet ween the observat i ons of iBraunl (Il995h an d simplified sunspot 
models were reported bv lFan et all (jl995h . lCallv et al.1 (|2QQ3h . and lCrouch et al.1 (|2005h 
using a forward modeling approach. Time-distance helioseismology and helioseismic 
holography aim at making images of the solar interior from maps of travel times or 
phase shifts under the traditional assumption that the Sun is weakly inhomogeneous 
in the horizontal directions. TD and HH h ave been used to infer wave speed varia- 
tions and flows in and ar o und sunspots (e. 



tions and hows m and around sunspots (e.g. 
Braun and Lindsevll2000l : iGizon et al1l200ol : 



Duvall Jr. et al. 1996; Jensen et al. 2001; 



M/ jKosovichev et al.1 l2000h IZhao et al.1 1200 ll : 
Couvidat etHI 2006) . Ring-diagram analysis has a coarser horizontal resolution and is 



used to study the s ubsurface structure of entire active regions (e.g. iBasu et al.ll2004l : 
Antia and Basu 20071; iBogart et al.ll2008T ). While these methods and their variants ap- 
pear to be quite robust, it has not been demonstrated that they are consistent. For 
instance ring-diagram analysis has not been directly compared to time-distance or 
holography in the case of an isolated sunspot. This paper reports on a joint study of 
the sunspot in NOAA Active Region 9787. 



2 Observations of NOAA Region 9787 

2.1 MDI/SOHO Observations 

NOAA Active Region 9787 was chosen from the MDI/SOHO data library because it 
hosts a large, round, isolated sunspot. A quick look at the data is given by Figure [T] 
The data consist of nine days of MDI full disk Dopplergrams for each minute from 20 
January to 28 January 2002. MDI also recorded the line-of-sight magnetic field every 
minute and intensity images every six hours. The images were remapped using Postel 
projection with a map scale of 0.12°. The centers of projection were chosen to track the 
motion of the sunspot (Carrington longitude cj) ~ 133° and latitude A = —8.3°). The re- 
mapping routine employs a cubic convolution interpolation. Missing data was linearly 
interpolated in time and a daily temporal mean was subtracted from each Dopplergram. 
Finally, we are left with one 512 x 512 x 1440 data cube of Doppler velocity data for each 
day. These data sets are made available on the European Helio- and Asteroseismology 
Network (HELAS) web site at http://www.mps.mpg.de/projects/seismo/NA4/, All 
authors were invited to analyse the same data, thereby eliminating discrepancies in 
the data reduction methods. 

Figure [T] shows a daily average of the MDI intensity continuum, magnetic field and 
Doppler velocity showing that there is little evolution of the sunspot during the period 
covered by the observations. The Dopplergrams show a ~ 2 km/s Evershed outflow in 
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Fig. 1 Daily averages of the SOHO/MDI Doppler velocity (top), intensity (middle), and line- 
of-sight magnetic field (bottom) of the sunspot in Active Region 9787 during January 20-28, 
2002. The Doppler velocity is in units of km/s, the magnetic field in units of kG. Each daily 
frame is a square with sides of length 200 Mm. 




Fig. 2 (Left) Normalised intensity profile of the sunspot averaged over time and azimuth. 
The vertical dashed lines indicate the umbral and penumbral boundaries. The boundary of 
the moat is given by the dotted line. (Right) Velocity of moving magnetic features (MMFs) 
averaged over time and azimuth as a function of distance from the center of the sunspot. 
The MMFs track the moat flow and are moving outward from the outer penumbra up to a 
radius of about 45 Mm (dashed line). The dotted indicate the boundaries of the umbra and 
the penumbra. 



the penumbra of the sunspot. The sunspot exhibits some amount of proper motion. 
Figure [2] shows the intensity profile of the sunspot averaged over nine days and over 
azimuthal angle, after correcting for the proper motion of the sunspot. The umbral 
and penumbral boundaries are at radii 9 Mm and 20 Mm respectively. 

The sunspot is surrounded by a region of horizontal outflow called the moat flow. 
In order to characterize the strength and extent of the moat, we measured the motion 
of the moving magnetic features (MMFs) from hourly averages of the magnetograms 
using a local correlation tracking method. The temporal and azimuthal averages of 
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the MMF velocity is plotted in Figure [2] as a function of distance from the center of 
the sunspot. The moat flow has a peak amplitude of 230 m/s and extends to about 
45 M m. The moat radius is about t wice the penumbral radius, which is a standard 
value ( Brickhouse and LabonteJll988h . 

The early development of Ac tive Region 9787 can be traced using the helioseis - 
mic technique of farside imaging (jBraun and Lindsev 2000; Li ndsev and Braunl hoOO). 
Figure [3] shows four maps of the full Sun created by K. Oslund and P.H. Scherrer (2006). 
These maps were taken from the SOI website at http : //soi . Stanford . edu/data/f ull_f arside/ . 
The overlaid grid represents Carrington solar coordinates where vertical lines of longi- 
tude are separated by 60° . Active Region 9787 is located close to latitude A = —8.3° and 
longitude <j) = 133° , as shown by the red circles in Figure [3] The active region is detected 
on the farside of the Sun in the top two panels of Figure [3] We then see the active region 
rotate past the East limb to the Earth side (third panel) and then within our observa- 
tion period (24th January 2002) in the bottom panel. Additional information is avail- 
able online at http://news-service.stanford.edu/pr/2006/pr-sun-031506.html. 



2.2 Oscillatory Power and Acoustic Halos 

There are two main properties of the acoustic power in and around active regions that 
are well documented. One is the power reduction in strong magnetic field regions, par- 
ticularly sunspots, and the second is the enhancement of power in t he higher frequencies 
(5-6 mHz) in the nearby photosphere, known as the acoustic halo (|Braun et alj fl992b: 
bonea et al. 2000). Here we quantify these properties for AR 9787. 

We calculate the temporal Fourier transform of the Doppler images for each day of 
observation. We divide this into 0.5 mHz bandwidths and calculate the power averaged 
over each of these frequency bandwidths. For all frequency bands the acoustic power 
suppression is greater than 80% in the umbra compared to the quiet Sun. We observe 
enhanced acoustic power at higher frequencies (5-6 mHz) in regions outside the sunspot 
and strong plage (Figure H}. 

Previous analysis of MDI acoustic power maps by lLadenkov et al.l (|2002h showed 
a modest excess of power around a sunspot in the higher frequencies, but which also 
appeared to be direct l y rela ted to the location of plage, rather than the sunspot itself. 
iHindman and Brownl (| 19981 ) also find that the high frequency velocity signal in an 
active region is higher (up to 60%) in pixels with moderate magnetic field strengths 
between 50 to 250 G. The fact that AR 9787 has an extended plage region offers 
an opportunity to analyze the plage far from the associated sun spot. In essence, this 
study bolsters the previous analysis of lHindman and Brownl {1 998) and Ladenkov et al.l 
(|2002h . 

When the sunspot is located close to the limb we find a significant enhancement of 
power in the umbra at high frequencies, to a level almost as high as in the quiet Sun. 
If this power is real and not an artifact of observing conditions at the limbs, then it 
could be due to magnetoacoustic waves with a component of motion perpendicular to 
the field lines in the umbra. 

From Figure [4] a region of particularly strong power close to the south-east side 
of the sunspot can be s een. This enh anced power is associated with the strong plage 
region just to the North. iBraunl (|l995T ) finds evidence of an acoustic deficit immediately 
outside a small (8 Mm) sunspot, extending out to 35 Mm and appearing to be well 
defined by the location of surrounding plage. The sunspot in AR 9787 shows little 
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Fig. 3 Maps of the full Sun created using helioseismic waves to infer the presence of magnetic 
activity on the farside of the Sun. Strong shifts (black/orange) in the phase indicate an active 
region, the light blue represents the quiet Sun. The horizontal lines are lines of constant 
latitude, the vertical lines are lines of constant Carrington longitude separated by 60°. The 
dates from top to bottom are 2002 June 15, 17, 19, and 24. Active Region 9787 is located close 
to latitude A = —8.3° and longitude <f> = 133° in all four maps, indicated by the red circles. 
(Courtesy of K. Oslund and P.H. Scherrer (2006)) 
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Fig. 4 Acoustic power averaged over all days and frequency bands 3-3.5 mHz (left) and 5.5- 
6 mHz (right). Overplotted contour is for B = 100 G, outlining regions of plage. The excess 
power outside the plage regions is clearly seen in the high-frequency maps (right). The power 
is normalised to unity in the quiet Sun (grey scale). 



evidence for a well denned, axisymmetric acoustic halo, leading to the suggestion that 
the enhanced acoustic power is associated with the strong surrounding plage regions, 
rather than the sunspot itself. 

The results of recent numerical work bv lHanasogel |2008) seem to reproduce the high 
frequency halo surrounding a small model sunspot. The halo appears b etween 5 and 
6 mHz, close to the acoustic cut-off. In agreement with the suggestion of iDonea et al.l 
(|2000h , the upwardly propagating fast mode waves may be reflected in magnetic regions 
due to the rapidly increasing Alfven velocity. These waves later re-emerge in the region 
surrounding the sunspot causing the observed enhanced power. 



2.3 Wave Absorption 

A useful analysis procedure for studying the interaction of p modes with sunspots is the 
decomposition of solar oscillations, observed in an annular region around the sunspot, 
into inward and outward propagating waves. Fourier- Hankel spectral decomposition 
has been used to identify p-mode absorption in sunspots and acti ve regions by com- 
paring the amplitudes of the ou tward and inward moving waves (|Braun et al.| [l988; 
iBogdan et al ] |l993l : lBraun1ll995h . 

Here we use a sunspot-centered spherical-polar coordinate system (0, 0) with the 
sunspot axis at = 0. The annular region is denned by the inner and outer circles at 
#min — 2.5° and #max = 11.25°, which correspond to distances between 30 and 137 Mm. 
The Doppler signal ip(0,<fi,t) in the annular region is decomposed into components of 
the form 

e i(m<t>+2n»t) ^ Am{ l^ )H W m + Bm{ l^ )H ^) m ^ , (1) 

where m is the azimuthal order, / is the harmonic degree, H^} and H^} are Hankel 
functions of the first and second kinds, t is time, v is temporal frequency, and A m 
and Bm are the complex amplitudes of the incoming and outgoing waves respectively. 
The range in / is between 70 and 1500. The boundaries of the annulus, 6 m [ n and #max, 



8 



c 4 




c 4 



200 400 600 800 1000 
Harmonic degree 



T — I — | — I — I — I — | — I — I — I — | — I — I — I — | — I — I — 



P 

r out 



■ 1.0 
J 0.9 
1 0.8 
1 0.7 
1 0.6 



"D 

10.5 i 



0.4 
0.3 
0.2 
0.1 



200 400 600 800 1000 
Harmonic degree 



Fig. 5 Power spectra of inward (left) and outward (right) propagating waves as a function 
of harmonic degree and frequency. The reduction in power of the outgoing modes is easily 
noticeable. 



were selected such as to resolve the low-order p mode ridges with a resolution in / of 
approximately 40. The numerical pro cedure needed to c ompute the wave amplitudes 
Am(l, v) and B m (l, v) is described bv lBraun et alj (|l98ct ). 

For each value of / we measure the mode amplitudes for the azimuthal order m — 0. 
The power spectra of the incoming and outgoing modes are displayed in Figure [5] The 
outgoing p-mode power appears to be significantly reduced compared to the incoming 
p-mode power. 

To measure this more quantitatively we determine an absorption coefficient for 
the / and p n ridges in a frequency band between 2.9 and 3.1 mHz. The absorption 
coefficient is defined as 



Jd/dz^ n (Pin-Pout) 
$dldvW n Pin 



(2) 



where W n is a window function that selects the n-th ridge, and P{ n (l, v) and P ut(l, v) 
are the power of the ingoing and outgoing waves. In this frequency band 2.9-3.1 mHz, 
the /, pi, p2, P3, and p^ ridges show absorption coefficients of 57%, 54%, 51%, 49%, 
and 50% respectively. These val ues are in agreement with those given in earlier studies 
(Bra un et al1ll988l : lBraun| [l995) and confirm that AR 9787 "absorbs" acoustic waves. 
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3 Travel Time Measurements 

3.1 Phase-Speed Filtering versus Ridge Filtering 

The principle of Helioseismic holography (HH) is to computat ionally regress the acous - 
tic amplitudes observed at the surface into the solar interior ([Lindsev and Braunll 19971 ) . 
To facilitate comparisons with results from time-distance analyses, we use surface- 
focused HH. In the "space-frequency" domain, i.e. where ip(r, v) denotes the temporal 
Fourier transform of the observed Doppler velocities, the regressions in surface- focused 
HH are computed from 

#± (r, v) = J dV G±(r, r', v) </>(r', v). (3) 

H+ and H— are the egression and ingression which represent estimates of the ampli- 
tudes propagating into and out of the focal point at position r on the surface and 
v is the temporal frequency. and G— are Green's functions that express how a 
monochromatic point disturbance at a position r on the surface propagates backward 
and forward in time, into the solar interior and back up to th e focus. They are computed 
using the eikonal approximation ( Lindsev and Braunlll997h . The correlations, 



C^{v) = {H^(v,v)r(v,v)) Av , (4) 

and 

C^(v) = ^(v,v)Hr*(v,v)) A „, (5) 

describe the egression and ingression control correlations respectively, which are di- 
rectly compa rable to center- annulu s correlations used in time-distance (TD) helioseis- 
mology (e.g. buvall Jr. et aljfliii : lBraunlll99/t ). The asterisk denotes complex con- 
jugation, and the brackets indicate an average over a chosen positive frequency range 
Av. 

Surface-focused HH can be used to study flows by dividing the pupil P, over which 
the ingressions and egressions are computed, into four quadrants (labeled N, S, E, 
and W), each spanning 90° and oriented in the North, South, East and West direc- 
tions respectively. We then compute the eight control correlations, (jN,S,E,W ^ Various 
combinations of these correlations are used to derive travel-time shifts due to the pres- 
ence of flows or wave speed perturbations. In general, we compute travel-time shifts 
from various sums or differences of correlations such that, if C denotes some linear 
combination of correlations, the travel-time shift is 

8r(r) = arg[C(r)]/2™ , (6) 

where vo is the central frequency of the bandpass Av. These represent travel-time 
shifts of the observed combination of waves relative to the travel times expected for 
the same ensemble of waves propagating in the solar model used to compute the Green's 
functions. 

We present measurements of the mean travel-time shift (^r m ean), which represents 
the shift computed from the sum of all eight correlations: 

<5r me a„(r) = arg[Cf (r) + Cf(r) + Cf (r) + C5(r) + 

(r) + (r) + Cf(r) + C^(r)]/2^ . (7) 
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We also present travel-time shifts sensitive to horizontal flows. For example, we define 
a "EW" travel-time asymmetry 

<5r E w(r, v) = (arg[C^(r, v) + (r, v )\ - arg[cf (r, v) + C™(r, u)])/2^ Q . (8) 

A similar travel-time shift can be measured for a North -South asymmetry. The sign 
of the travel-time perturbations is such that a perturbation in the background wave 
speed resulting in a faster propagation time will lead to a negative value of the mean 
travel-time shift (^r m ean) and a horizontal flow directed from East to West produces 
a negative value of the EW travel-time asymmetry (<5tew)- 

Starting with the tracked, Postel-projected datacube described earlier fSection !2.ip 
we perform the following steps: 1) a temporal detrending by subtraction of a linear fit 
to each pixel signal in time, 2) removal of poor quality pixels, identified by a five- 
sigma deviation of any pixel from the linear trend, 3) Fourier transfor m of the data in 
time , 4) a correction for the amplitude suppression in magnetic regions ([Raiaguru et al.l 
2006), 5) spatial Fourier transform of the data and multiplication by a chosen filter, 
6) extraction of the desired frequency bandpass, 7) computation of Green's functions 
over the appropriate pupils, 8) computation of ingression and egression amplitudes by 
a 3D convolution of the data with the Green's functions, and 9) computation of the 
travel-time shift maps by Equations Q - J5} . 

In step 5 we have used two general kinds of filters: phase-speed filters and ridge 
filters. While not as commonly use d as phase-speed filters, ridge filters have been used 
previously for /-mode studies (e.g. iDuvall Jr. and Gizo n 2000: lGizon and Birchll2002l: 



Jackiewicz et ajj 2007; I Jackiewicz et al.ll2007h and recently for p modes ([jackiewicz et alj 



2008; iBraun and Birchl l2008h . The phase-speed filters used here are the s ame set of 
11 filters (here a fter d enoted "TD1-TD11") listed bv iGizon and Birchl (|2005h and 



Couvidat et~ai1 (|2006h and commonly used in time-distance analyses. The inner and 



outer radii of the corresponding pupil quadrants are chosen so that acoustic rays at a 
frequency of v — 3.5 mHz propagating from the focus to the edges of the pupil have 
phase speeds (denoted by w) which span the full width at half maximum (FWHM) of 
the squared filter. The EW pupil quadrants for the 11 filters are shown in the bottom 
set of panels in Figure |6] All of the pha se speed filters used also remove the contribution 
of the / mode (jBraun and Birchl l2008h . 

Figure [6] shows mean and EW travel-time shifts using a frequency bandpass (step 
6) between 2.5-5.5 mHz. As is well known, there are distinct patterns of travel-time 
shifts associated with the use of phase-speed filters. In particular, the mean travel- 
time shift in the sunspot is positive for the smallest values of phase speed (and mean 
pupil diameter; e.g. TD1-TD3) and then switches sign for larger phase speeds (TD4 
and beyond). The EW travel-time shifts also undergo a similar change of sign. At 
the smallest (largest) phase speeds, the EW travel-time differences are consistent with 
inflow- like (outflow- like) perturbations centered on the sunspot. Remarkably, the switch 
in sign for both the mean shift and the EW differences occurs between filters TD3 and 
TD4. 

Figure [7] shows mean and EW travel-time shifts using the same frequency bandpass 
(2.5-5.5 mHz) as shown in Figure[6l but obtained with ridge filters isolating the p\ — p4 
ridges. In contrast with the results obtained using phase-speed filters, ridge filters show 
values of £r m ean which are always negative within the sunspot, while £tew is consistent 
with an outflow-like perturbation . This is similar to resul ts from time-distance analyses 
as well as previous HH analyses ( Braun and Birchl l2008h . 
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Fig. 6 Maps of travel-time shifts Sr me an (top panels) and 5tew (bottom panels) using phase 
speed filters over a frequency bandpass of 2.5-5.5 mHz for AR 9787 observed over a 24 hour 
period on 2002 January 24. The labels TD1 through TD11 indicate the phase-speed filter used. 
Sizes of the East and West pupil quadrants used to measure <5tew are shown in the bottom set 
of panels. The map in the lowest-left position of the top set of panels shows a MDI continuum 
intensity image while the map in the same position in the bottom set shows a line-of-sight 
magnetogram of the sunspot in AR 9787. The portion of the region shown here extends 219 
Mm on each side. For the purpose of this Figure (and the others in this section), some spatial 
smoothing has been applied to the travel-time maps. 
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Fig. 7 Maps of travel-time shifts 5r m e&n (top panels) and Sr^w (bottom panels) using ridge 
niters over a frequency bandpass of 2.5-5.5 mHz for AR 9787 on 2002 January 24. The columns 
indicate different radial orders as indicated. 

Motivated by recent s tudies which show frequency variations of travel-time shifts 

obser ved in active regions (|Braun and Bir ch 2006; Couvidat and Raiagurull2007l ; lBraun and Birchl 

2008), we also employed narrow frequency bandpasses to both the phase-speed and 

ridge filters. Figure [8] shows mean and EW travel-time shifts using phase-speed niters 

TD1-TD4 in conjunction with 1-mHz wide frequency filters centered at 2, 3, 4, and 

5 mHz. Positive mean travel-time shifts, and EW travel- time differences consistent 

with inflows, are observed primarily in frequency bandwidths that are centered below 

the pi ridge, shown by the solid line in Figures [8] There is one instance (filter TD3 at 

5mHz, which lies immediately below the p2 ridge) which also produces a positive mean 

shift and inflow- like signature. All other filters (including TD5-TD11 not shown) show 

nega tive mean travel - time s hifts and outflow-like signatures. 

iBraun and Birchl (|2008h find that a condition for producing positive travel-time 
shifts such as in Figure [8] appears to be a disproportionate contribution to the cor- 
relations of wave power from the low- frequency w ing of the p\ ridge relative to the 
high-frequency wing. We note that recent work bv lMoradi et al.l (|2008h is relevant to 
these issues. 

3.2 Ridge and Off-Ridge Filtering 

Follow ing IBraun and Birchl (|2006h , lBraun and Birchl (|20Q8h and lThompson and Zharkovl 
(2008), we study the sensitivity to filtering of mean travel-time perturbations measured 
in the vicinity of isolated sunspot AR 9787 relative to the surrounding quiet Sun, using 
a centre-to-annulus geometry and a skip distance of 11.7 Mm (Figure For each row a 
bandpass filter was used to select data within a 1 mHz frequency band with a 0.1 mHz 
Gaussian roll-off, centred (from bottom to top) at 2.5, 3.0, 3.5, 4.0, 4.5 and 5.0 mHz. 
These were combined for each column (left to right) with filters selecting the data from 
in between / and p\ ridges, from the p\ ridge, in between the p\ and p2 ridges, and 
the p2 ridge. The filters were constructed as follows: at constant frequency we apply 
a filter that takes the value of unity at the horizontal wavenumber corresponding to 
either a particular ridge, e.g. p\ (a "ridge filter"), or a mid-point between the adjacent 
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Fig. 8 Maps of travel-time shifts Sr me an (top panels) and 5tew (bottom panels) using phase 
speed filters and 1 mHz-wide frequency bandpasses for AR 9787 observed on 2002 January 24. 
The columns of maps labeled TD1 through TD4 indicate the phase-speed filter used, while 
the rows indicate the frequency bandpass. The solid jagged line running diagonally through 
the panels connects the location of the pi ridge in the u-w domain for each filter, with the 
centers of the maps assigned to values of frequency and phase speed as indicated on the left 
and bottom edges of the plot. The dashed and dotted lines indicate the locations of the p2 
and p3 ridges respectively. The map in the lowest-left position of the top set of panels shows a 
MDI continuum intensity image while the map in the same position in the bottom set shows 
a line-of-sight magnetogram. 
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ridges, e.g. p\ — p2 (an "off-ridge filter"). On either side of this centre line the filter has 
a Gaussian roll-off with half width at half maximum (HWHM) equal to 0.32 times the 
distance to the neighbouring ridge on that side for the ridge filter and with HWHM 
equal to 0.63 times the distance to the adjacent ridge in the case of the off-ridge filter. 
No phase-speed filter was applied. 

The cross-correlation functions were estimated using center-to-annulus geometry 
with annuli taken to b e one data pixel wide. T ravel-time perturbations were measured 
using the definition of iGizon and Birch 

In agreement with Braun and Birchl (|2006h and iBraun and Birchl (|2008h , as illus- 
trated in Figure [9] we observe a positive travel-time perturbation in the region be- 
neath the pi ridge, but we also find such a signal between the p\ and p2 ridges, and 
find that on the p\ ridge the positive perturbation is absent. Our tentative conclusion 
is that the positive travel-time perturbation signal arises only in the regions between 
the p-mode ridges; and that the travel-time perturbations associated with the data 
on the ridges themselves are all consistently negative. Similar, though noisier, results 
were obtained for the Gabor- wavelet fitting travel-time definition. The data appear to 
suggest that measured travel-time perturbations are very sensitive to the part of the 
wave-propagation diagram selected during the filtering stage. 



4 Moat Flow Inversion (Ridge Filters) 

Here we invert travel times to obtain the flows around the sunspot in AR 9787. We 
restrict ourselves to ridge filtering and TD travel times. We then compare the inferred 
flows with the velocities of the moving magnetic features (MMFs) in the moat (Sec- 
tion l2Tl). 

The details of the travel-time inversion can be found in ljackiewicz et all (|2QQ8h . In 
summary, we measure center- to- quadrant travel-time differences using the method of 
IGizon and Birchl (2004), after correcting for reduced power in magnetic r egions. Using 



the s ame definition of travel times, Born sensitivity kernels are computed ( Birch and Gizonl 



l2007h . We input the travel times and the sensitivity functions, as well as the covariance 
of the travel times, into a three-dimensional subtractive optimally localized averages 
(SOLA) inversion procedure to infer the vector flows at several depths. The procedure 
also provides good estimations of the resolution and the noise levels, which are impor- 
tant for any interpretation. We note that neither the modeling nor the inversion takes 
into account the magnetic field. 

We have inverted the middle seven days of data from the nine day set. We obtain 
flow maps for several depths, extending down to about 5 Mm beneath the surface. 
In Figure [10] we compare the inversion results near the surface (left column) to the 
motion of the moving magnetic features (MMFs, right column). For both sets of maps, 
the flows are averaged over 7 days. Furthermore, they are approximately of the same 
horizontal resolution (~ 6 Mm). For this particular figure, we choose to study the 
inferred time-distance flows taken at a depth as near to the surface as we can achieve, 
about 1 Mm below. 

The bottom row of Figure [10] compares the radial velocities derived from TD he- 
lioseismology and MMF tracking. We see for each case quite clearly a strong outflow 
extending beyond the penumbra (20 Mm) of several hundred m/s, known as the moat 
flow. The overall features of the flows from both methods are quite similar, even the 
slight 'knob' on the northeast quadrant of the sunspot moat. A look in the quiet Sun 
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Fig. 9 Travel-time perturbations 5r m ean for isolated sunspot AR 9787 obtained for various 
filtering schemes. The colorbar is in units of seconds. The skip distance is equal to 11.64 
Mm. From the left, the columns respectively have niters applied as follows: pass-filter centred 
between the / and pi ridges; pass-filter centred on the pi ridge; pass-filter centred between the 
pi and p2 ridges; pass-filter centred on the p2 ridge. From the bottom, the rows have bandpass 
filters centred on 2.5, 3.0, 3.5, 4.0, 4.5, and 5.0 mHz respectively. More details are given in the 
text. 



reveals other similarities. The correlation coefficient between the two maps is about 
0.65. The magnitude of the TD surface flows is about 20% less than that revealed by 
MMF tracking. This can be due to many factors, such as the implied depth at which we 
are comparing not being equal, inaccurate travel-time sensitivity kernels, or the mag- 
netic field affecting the inversion results through the travel times, among others. The 
estimated noise in the time-distance maps is about 5 m/s. These results are consistent 
with a previous /-mode TD study of sunspot moat flows (|Gizon et al.ll2000h . 

Since we do not know the exact depth where the small magnetic features are rooted, 
it is useful to compare their motion with inverted flows at several depths. Since the 
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Fig. 10 Comparison of near-surface flows around the sunspot from time-distance inversions 
and MMF tracking. The left column shows flows obtained from inversions at a depth of 1 Mm 
beneath the surface and averaged over 7 days. The right column are the flows obtained from 
MMF tracking averaged over the same 7 days. Each set of flows was obtained with approx- 
imately the same resolution/smoothing per pixel. The top row is the x component of the 
velocity, and the middle row shows the y component of the velocity. The bottom row compares 
the radial velocity (from the center of the sunpot) of the two measurements. An outward moat 
flow is seen for each case. 
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Fig. 11 Azimuthally-averaged radial flows from the sunspot center at different depths over 
7 days, obtained from a time-distance inversion. The radial velocity from the MMF tracking 
is shown by the dashed line. In this plot only the region from the edge of the penumbra 
outward is shown. The results within the sunspot (distance less than 20 Mm) cannot be 
reliably interpreted. 

sunspot is very nearly circular, we may first average the radial flows azimuthally about 
the sunspot center. In Figure [TT] we plot the azimuthally-averaged radial flows for 
several depths versus the distance from the center of the sunspot. Also shown for 
comparison is the averaged MMF velocity. The moat flow, which extends to about 
45 Mm, is seen at all depths. The inferred flows get stronger with depth and begins to 
level off in strength at about 4.5 Mm below the surface. The MMF velocity is consistent 
with these flows at depths of between 1 and 2.6 Mm beneath the surface. The shape 
of the moat flow is similar for all cases. 

5 Sound Speed versus Wave Speed 

The sound speed, c, is the speed of sound. In this paper we make an important dis- 
tinction between c and the speed at which waves are inferred to have travelled under 
the assumptions of the helioseismic inversions. Different inversions will include differ- 
ent physics. A common example is the assumption that a travel-time perturbation 
associated with a sunspot can be modelled purely by an equivalent small-amplitude 
sound-speed perturbation. The sound-speed inferred under such an assumption is at 
best a local wave speed, and need not even reflect the sign of the sound speed pertur- 
bation. For this reason, we use the notation c w to denote the inferred wave speed. 
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I Lin et all (|2008h studied the exact meaning of c w in the case of ring-diagram anal- 
ysis. They find that their inversions return perturbations in = Aftot/p, where 
Ptot — P + -Rnag is the sum of the gas and magnetic pressures. The inferred wave 
speed c w has two components: the sound speed, c, and a magnetic component. 

What we have said about sound speed versus wave speed, also applies to any chosen 
physical quantity versus inferred quantity. For example, the next section presents an 
inversion for wave slowness, s w = l/c w , using traveltime sensitivity kernels for inverse 
sound speed s = 1/c. 



6 Wave-Speed Inversion (Phase-Speed Filters) 

The data was reduced using standard time-distance helioseismology techniques. First, 
th e data was preprocessed by applying an amplitude modulation correction as described 
bv lRaiaguru et al.1 (200(5). followed by applying a high-pass filter at 1.7 mHz in order to 



remove the supergranulation and a low-pass filter at 5.1 mHz to remove signal above the 
acoustic cut-off frequency. We then apply a Gaussian phase speed filter to select waves 
with horizontal phase speed W{ near the value corresponding to skip distance A{ given 
by ray theory. We define 12 such filters Fi corresponding to different distances Aj. These 
phase-speed filters are similar to the set of filters used in Section 13.11 and commonly 
used in time-distance analyses. Each filter is applied by pointwise multiplication of the 
Fourier transform ip(k x , k y ,uj) of the observed velocity data. 

We use a centre-to-annulus geometry to compute the cross- covariance C(r, Aj,£), 
where r is the center of the annulus and t is the time lag. The annulus width is 4.5 Mm. 
A reference cross- covariance function, C ref , is obtained by spatial averaging C over 
a quiet Sun area. Wave travel times are then ex tracted by fitting a Gabor w avelet 
to the positive- and negative-time branches of C ( Kosovichev and Duval]||l997h . The 



wavelet has five parameters: the central frequency, the width and amplitude of the 
envelope, and the group and phase travel times. We denote by r+ and r_ the measured 
phase travel times for the positive- and negative-time branches of C respectively. The 
reference travel times for the quiet Sun are similarly defined using C ref . The phase 
travel time perturbations, 5t+ and 6t—, are defined as the difference between the 
measured and reference travel times. As we are interested in wave-speed perturbations 
only, we consider mean travel-time perturbations, ^r m ean = (£t+ + St—)/2. 

For the forward p roblem we use sensitivity kernels estimated using the first-order 
Rytov approximation (| Jensen and Piipersll2003r ). These kernels, K s , relate mean travel- 
time perturbations, ^r m ean, to inverse sound speed perturbations, 5 s — 5 (1/c), of a 
quiet-Sun model. In the sunspot region, we have 

Sr m ean(r,A)= f d 2 r' f dz K s (r — r', z; A{) 6s w (r , z), (9) 
Js J-d 

where S is the area of the region, d is its depth. The quantity 5s w — Sw S IS the 
equivalent change in the local wave slowness caused by the sunspot. 

We invert for N = 14 layers in depth located at [z±, • • • zjy] = [0.36, 1.2, 2.1, 3.3, 
4.7, 6.4, 8.6, 11.2, 14 .3, 17.8, 21.8 26.3 31.4 37.01 Mm. We use a multi-channel decon- 
volution algorithm (Ijensen et al.l Il998l . l200lh enhanced by the addition of horizontal 



regularization ([Couvidat et al. l2006h . The above equation is Fourier transformed with 



respect to two-dimensional position r. For each wavevector k, we define di = 5r(k, 



19 



u 1 ° 

CD 



15 




-100 -50 50 100 
Mm 

Fig. 12 Relative wave speed perturbation, 5c w /c, obtained for AR 9787 using time-distance 
helioseismology and phase speed filtering. 



and rrij = Ss(k, zj) 



and the corresponding vector d, matrix G, 



and vector m. Then for each k we solve for the vector m that minimizes 

||(d-Gm)|| 2 + e||Lm|| 2 , 



(10) 



where L is a regularization operator and e(k) is a positive regularization parameter. 
In this work we apply more regularization at larger depths, to which travel times are 
less sensitive, by setting L = diag(ci, C2, . . . cjy), where cj = c(zj) is the sound speed 
in the j-th layer of the reference model. We regularise small horizontal scales by taking 
e(k) = 2 x 10 3 (1 + |k| 2 ) 100 . 

Figure [12] shows the result of the inversion, expressed in terms of the relative 
wave-speed perturbation Sc w /c. We see a two-layer structure: a region of decreased 
wave speed (down by —13%) situated directly underneath the surface and a region of 
increased wave speed (up to 9%) starting from a depth of approximately 3 Mm. This is 
consistent wi th other time-dist a nce in versions of travel times using phase-speed niters, 
e.g., those of iKosovichev et all ([2000h . 



7 Ring-Diagram Analysis 

We have used ring diagrams to analyze the mean structure of the region containing 
A R 9787 compared w ith quiet-Sun structure. To do so, we use the techniques described 
in iBasu et al.l (|20Q4h . We invert the differences in the ring-diagram fit parameters 
between the spectra of the active region and those of suitable selected quiet- Sun regions. 
In this case, two quiet-Sun regions were chosen at the same latitude as that of the 
active region (—7°) and at Carrington longitudes 170° and 75° (the active region is 
at longitude 130°). Each region is independently tracked in a time interval of 5.7 days 
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Fig. 13 SOLA inversions for the depth dependence of the relative difference in wave speed, 
5c w /c, between the region around AR 9787 and each of two comparison quiet regions at the 
same latitude but different longitudes. The inversions are based on fits to power spectra for 5.7 
days of tracked data for each region. The differences are in the sense Sc w = c W)ac ti ve — Cw, quiet- 
The active region is seen to have a positive wave-speed anomaly very near the surface relative 
to the quiet regions. 



centered on its central meridian crossing, so any geometrical differences in the spectra 
due to foreshortening or geometric image distortion are very nearly canceled out. The 
regions chosen for analysis are 15° in diameter, so the results apply to a spatial mean 
over these areas (with an unknown weighting function). The average results from the 
two comparison regions is shown in Figure[l3] It is clear that there is a region of negative 
wave-speed anomaly under the active region between 3 Mm and 8 Mm in depth, with 
a turnover to positive wave-speed anomalies both above and below this region, and 
yet another turnover to negative anomalies at depths gre ater than about 17 Mm. This 
behaviour is typical of that seen for other active regions (|Bogart et al.ll2008r h although 
the changes at the surface and deeper than 17 Mm are unusually pronounced in this 
case. 
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Fig. 14 SOLA inversions for the depth dependence of the mean zonal (u x ) and meridional 
(u y ) flows over each of the three regions analyzed with ring-diagram analysis. The zonal rates 
are relative to the tracking velocity, which was 54 m/s less than the Carrington rate at this 
latitude. 



We can also infer the mean flow speeds at depth for the active regions and the 
comparison quiet regions directly from the fitted parameters to the ring-diagram spec- 
tra. These are shown in Figure ITU There is no evident anomalous zonal flow through 
the active region; indeed, the zonal flow structure is remarkably similar to that of the 
preceding comparison region at longitude 170°. There does appear to be an anomaly 
in the structure of the meridional flow, however, with a substantial shear at depths 
greater than 7 Mm, the flow being poleward near the surface (the region is in the 
southern hemisphere) and equatorward at greater depths. It is especially marked if the 
mean meridional velocity at the active region's latitude is negative at depth, as the 
two comparison regions suggest, but this needs to be verified by averaging over more 
longitudes. 
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Fig. 15 Comparison between the observed cross-covariance and a numerical simulation. The 
circle indicates the location of the sunspot in AR 9787. (Top half) /-mode cross-covariance 
between the MDI Doppler signal averaged over a line at x = —40 Mm and the Doppler signal 
at each point. The correlation time-lag is 130 min, large enough for wave packets to traverse 
the sunspot. The cross-covariance is averaged over 9 days and uses the assumed azimuthal 
symmetry of the sunspot to reduce noise. (Bottom half) SLiM numerical simulation of an 
/-mode wave packet propagatin g in the -\-x direction through a model of AR 9787 with a peak 
magnetic field of 3 kG. (Trom I Cameron et al.l (|2008l ^ 

8 Numerical Forward Modeling 

In various circumstances it has been shown that the cross-covariance is closely related 
to the Green's function. This allows us to characterize the interaction of arbitrary 
wavepackets with the sunspot from the MDI observations. The sunspot discussed in this 
paper, being observed over nin e days and almos t axis ymmetric, is ideally suited to such 
a study. In preliminary work, I Cameron et al.l (|20Q8h considered the cross-covariance 
between the Doppler signal averaged along a great circle 40 Mm from the centre of the 
sunspot and the Doppler signal at each point in a region surrounding the sunspot. The 
data had been /-mode ridge filtered. 

A numerical simulation was then performed of the propagation an /-mode plane 
wave packe t beginning 40 Mm from a model s unspot. The background atmosphere is 
Model S of Christensen-Dalsgaard et al. I CHI), 

stabilized with resp ect to convection. 

The su nspot model used was a simple self-similar model in the vein of lSchliiter and Temesvarvl 
(|l958lb . The half- width of the vertical magnetic field at the surface was taken to be 
10 Mm (as in A R 9787) and field stre ngths of 2000, 2500 and 3000 G were considered. 
The SLiM code ([Cameron et alj 12007) was used to perform the simulations. 

The top half of Figure[T5lshows the observed cross-correlation with time-lag 130 min. 
The bottom half of the figure shows the results of the simulation at time 130 min. The 
match between the observations and simulation is quite good, in this case for a sunspot 
with a peak field strength of 3000 G. The match was not as good for peak field strengths 
of 2000 G and 2500 G. This then places a helioseismic constraint on the magnetic field of 
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Fig. 16 Structure of the 3 kG sunspot model used in the simulation of Figure El The 
thick solid curve shows the relative change in gas pressure, P, with respect to the quiet Sun, 
measured along the sunspot axis. The thick dashed curve shows the relative change in the 
squared sound speed, c 2 , which is also an estimate of the relative change in temperature. The 
vertical dashed line indicates the depth at which the sound speed and the Alfven velocity are 
equal. The vertical dotted line indicates the depth at which the sound speed is equal to its 
zero-depth quiet-Sun value. 



the spot. Whereas this constraint makes sense, it cannot be assessed directly using MDI 
magnet ograms, which are not reliable in sunspot umbrae. The partial absorption of the 
waves (reduced cross- covariance amplitude) was explained in terms of mode conversion 
into slow magnet o-acoustic waves that propagate dow n the sunspot, as pred icted by, 
e.g., [CallyJ (2000). Full details of this work are given in I Cameron et ahl (12008). 



For the 3 kG sunspot model, Figure [TBI shows the relative change in the gas pressure 
and the square of the sound speed along the sunspot axis with respect to the quiet Sun 
model. Both perturbations become very small (less than 1%) deeper than a depth of 
about 4 Mm. The temperature, closely related to the squared sound speed, is reduced 
at all depths within the sunspot. At the surface the relative decrease in temperature is 
around 18%. The reduction in gas pressure is larger with a 36% decrease. The vertical 
dashed line indicates that the sound speed is equal to the Alfven speed (c = a) at a 
depth of approximately 580 km. The c— a level is where mode conversion is expected 
to occur. We also plot the depth at which the sound speed is equal to the surface 
quiet-Sun sound speed. This gives a rough indication of the Wilson depression, in this 
case a rather low 170 km. 
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9 Discussion: The Elusive Structure of Sunspots 

9.1 Problematic Travel Times 

Here we use Born approximation based forward modeling to test if the flows estimated 
from inversions of the travel-times obtained using ridge filters (Section [4} are consistent 
with the travel-times measured using phase-speed filters (Section [3]). 

To carry out the forwa rd modeling we employ sensitivity functions, kernels, com- 
puted using the method of iBirch and Gizonl (J2007|). The calculations account for both 
the phase-speed filters and pupil sizes used for the measurements described in Section[3] 
The resulting kernels, K v , relate three-dimensional steady flows, v, to predictions for 
the EW travel-times differences, 



where r and r' are two-dimensional position vectors and z is depth. Both the kernel 
functions and the flow are vector- valued functions of horizontal position and depth. The 
travel-time differences are functions of horizontal position and also the pupil size A 
(notice that for each pupil size there is a corresponding phase-speed filter, see Section[3j). 

We assume v to be given by the flow field inferred from the inversions of the ridge- 
filtered travel-time differences shown in Section [U We neglect the effects of vertical 
flows as inversions for depth-dependent vertical flows have not yet been carried out. 

Figure [T71 compares measurements and forward models of EW travel-time differ- 
ences for the cases of phase-speed filters TD1 and TD4. For the case of the filter TD1, 
the forward model is in qualitative agreement with the measurements in the quiet Sun. 
However, the forward model predicts a signature of the moat flow with opposite sign 
to that seen in the observations. This shows that the travel-time differences measured 
using phase-speed filter TD1 and those measured using ridge- filtered travel-times do 
not yield a consistent picture of the moat flow. Figure \T7\ also shows the case of mea- 
surements and forward modeling for phase-speed filter TD4. In this case, and for other 
filters with large phase speeds, there is qualitative agreement even in the moat. 

One possible reason for the disagreements could be that the travel-time sensitivity 
kernels rely on an imperfect model of the power spectrum of solar oscillations. For 
example, the model zero-order power spectrum does not include background noise and 
mode linewidths may not be accurate enough. 

We also emphasize that a number of assumptions have been made in carrying 
out the forward model i ng shown in Fig ure 1 171 It is known (e.g. iGizon and Birchlkocd 
IParchevskv et all 20081 ; Irlanasoge et al.l2007h that the reduction in the wave generation 
rate in sunspots can, in general, produce apparent travel-time differences. For the case 
of p hase-speed filtered tra vel-time differences this effect has a magnitude of up to 
10 s (jHanasoge et al.ll2007f ). The magnitude of this effect is not known for travel-times 
measured using ridge filter s. Similarly, wave damping in sunspo ts can also produce 
travel-time differences (e.g. IWoodardlll997l ; IGizon and Birchi r2002 N ). The magnitude of 
this effect has not been carefully estimated for realistic models of wave absorption in 
sunspots. In addition , radiative transfer effects can cause phase shifts in sunpots (e.g. 




(ii) 



iRaiaguru et al.ll2007t ). 
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Fig. 17 Measured EW travel-time differences for phase-speed niters TD1 (panel a), and 
TD4 (panel c), and the corresponding modeled travel times (panels b and d) for 25 January 
2002. The travel-time measurements are described in Section [3] and have been smoothed with 
a Gaussian filter with FWHM of 6.7 pixels. The color bars have units of seconds. The gray 
scales are shown in units of seconds and have been truncated to make the details more visible. 
The sunspot is centered at roughly (x,y) = (0,0). Notice that the forward model is able to 
reproduce the moat flow for the filter TD4, but not for the case of TD1. 



9.2 Conflicting Wave- Speed Profiles 

Here we compare the wave-speed inversions from ring analysis (Section [7J) and from 
time-distance helioseismology with phase-speed filters (Section [6]). To make this com- 
parison possible, we average the TD wave-speed inversion over a disk of 15° diameter 
centered on the sunspot, which represents the area used for ring analysis. The two 
wave-speed profiles are plotted as a function of depth in Figure [18] Clearly, they do 
not match. 

How can we explain such a strong disagreement? As already mentioned in Sections[3] 
and 19.11 the details of the measurement procedures are important for the interpretation 
of the helioseismic observations; they may not have been fully taken into account in 
one or possibly both inversions. Although we have not done a TD inversion for wave 
speed using ridge filters, it is likely that it would give a different answer than the TD 
inversion using phase-speed filters, thus adding a third curve to Figure [T8l 

We also note that both inversions suppose that first-order perturbation theory is 
valid to describe the effect of sunspots on waves. Unlike the flow perturbation, how- 
ever, the perturbations in pressure and density introduced by the sunspot are not small 
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Fig. 18 Comparison of two different helioseismic methods used to infer wave speed per- 
turbations below AR 9787 (Sc w /c). The red curve shows the averaged ring-diagram results 
from Figure fT3l The solid blue curve shows the time-distance result (phase-speed filters) from 
Figure 1121 after averaging over the area used for ring analysis. Although they are meant to 
represent the same quantity, these two curves are noticeably different. 

with respect to the quiet-Sun background. Thus the concept of linear inversions is not 
necessarily correct for sunspots and regions of strong magnetic field. In addition, it is 
perhaps too naive to model the combined effects of the magnetic field in terms of an 
equivalent sound-speed perturbation. Ring-diagram inversions do include a contribu- 
tion from changes in the first adiabatic exponent, but the direct effect of the magnetic 
field through the Lorentz force is not fully accounted for in either inversion. We note 
that the ring-diagram inversions include a treatment of near-surface effects which is 
different than in the TD inversions. 



10 Conclusion 

We have studied the sunspot in AR 9787 with several methods of local helioseismology. 
We have characterized the acoustic wave field near the sunspot and the surrounding 
plage, measured acoustic absorption by the sunspot, and showed maps of the signature 
of AR 9787 on the farside of the Sun. We have shown that the sunspot leaves a strong 
signature in the observed wave field, as evidenced by strong perturbations in travel 
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times and frequency shifts. The interpretation of the observations, however, is difficult 
and we have not been able to draw an unequivocal conclusion about the subsurface 
structure and dynamics of the sunspot. We have shown that one complication is the 
extreme sensitivity of helioseismic measurements to the choice of data analysis proce- 
dure, such as filtering in frequency- wavenumber space. In addition, our understanding 
of the effects of strong magnetic fields on solar oscillations is still incomplete. 

On the positive side, we note that the seismically determined moat flow (TD and 
ridge filters) appears to be consistent with the motion of the MMFs in magnetograms. 
It is also clear that numerical simulations of wave propagation through model sunspots 
promise to provide invaluable help in interpreting the observations. 
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